Squeezing the Collective Spin of a Dilute Atomic Ensemble by Cavity Feedback 
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We propose and analyze a simple method to squeeze dynamically and unconditionally the col- 
lective spin of a dilute atomic ensemble by interaction with a driven mode of an optical resonator, 
as recently demonstrated [I. D. L., M. H. S., and V. V., Phys. Rev. Lett. 104, 073602 (2010)]. 
We show that substantial squeezing can be achieved in the regime of strong collective ensemble- 
resonator coupling. The squeezing is ultimately limited either by photon emission into free space or 
by the curvature of the Bloch sphere. We derive both limits and show where each prevails. 
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While techniques for preparing single-particle spin 
states in atomic ensembles (coherent spin states are 
well established, the manipulation of arbitrary quantum 
mechanical many-body states remains far out of reach, 
and only a tiny fraction of a many-spin Hilbert space 
is experimentally accessible to date. The preparation of 
even weakly quantum-correlated (entangled) states re- 
tires a controllable interaction between the particles 
that induces system evolution that is fast compared 
to system decoherence. Only a few (pseudo-)spin sys- 
tems offer such favorable interaction-to-decoherence ra- 
tios, notably trapped ions, where entangled states of up 
to eight particles have been prepared S 4| , and colliding 
atoms in a Bose- Einstein condensate [5|- In those sys- 
tems, squeezed spin states [U, S 0] — m which the quan- 
tum noise is redistributed so that one noise component is 
smaller than possible for unentangled states — have been 
generated [8|, |9j- 

In a dilute atomic ensemble, a common interaction of 
the atoms with a light field can replace direct spin-spin 
interactions 

EMU- Entanglement between the spin and 
light degrees of freedom allows a measurement performed 
on the light to reduce spin noise [ill [l3j and produce spin 
squeezing conditionally, as recently demonstrated by two 
groups using trapped atom s fbH 151. 

Pioneering proposals [16- 18] have shown that repeated 
light-ensemble interaction can produce spin squeezing dy- 
namically and deterministically, without requiring mea- 
surement of the light field. Such processes can be viewed 
as quantum coherent feedback [l9| : the ensemble spin im- 
prints its quantum fluctuations on the light, which acts 
back on the spin state to reduce those fluctuations. 

Here we propose a simple and robust method for 
squeezing the spin of an atomic ensemble inside an optical 
cavity by coherent feedback. We show that a laser tuned 
to the slope of the cavity resonance and off-resonant 
from the atomic transition induces a one-axis twisting 
@; 13 °f the spin state space (Bloch sphere; Fig. Id). 
For an initial coherent spin state (CSS) Q orthogonal 
to the twisting axis, this reduces the quantum noise in 
an oblique spin component. While photon shot noise in 
the incident light and photon scattering into free space 
counteract the noise reduction, substantial squeezing is 
possible provided the collective cooperativity (resonant 
optical depth) is large, even if single atoms are weakly 
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FIG. 1: (Color online) Scheme for spin squeezing by cavity 
feedback. An ensemble of atoms is uniformly coupled to a 
laser held in an optical resonator (a). A laser tuned to the 
slope of a cavity resonance (b), at equal and opposite detun- 
ings from transitions |f) — > |e) and |4-) —¥ \e) (c), shears a 
coherent spin state prepared along the x axis into a squeezed 
spin state (d), illustrated by a tomographic probability dis- 
tribution [2y| on the Bloch sphere. 



coupled to the cavity. 

For an ensemble of two-level (spin-^) atoms, described 
by a collective spin S, the squeezing is produced by an 
ensemble-light interaction Hamiltonian of the form c'cS z 
that represents the differential energy shift between the 
two atomic states due to the intracavity light with pho- 
ton number c^c. This photon number depends on the 
population difference 2S Z between the atomic states be- 
cause the precise tuning of the resonator mode relative to 
the driving laser depends on the atoms' state-dependent 
index of refraction. In particular, if c^c depends linearly 
on S z , the differential energy shift causes a precession of 
the spin vector about the z axis at a rate proportional 
to S z , similar to the dynamics of a Hamiltonian oc S z 
[2|. The initially circular uncertainty region of a CSS is 
then sheared into an ellipse (Fig. [TJl) that is narrower 
in one direction than the original CSS, corresponding to 
spin squeezing. The spin correlations between different 
atoms arise from the fact that the phase between the two 
states in any individual atom now depends on the state 
population difference 2S Z of the entire ensemble. 

We consider an ensemble of N identical three-level 
atoms whose ground states |f) , \i) (e.g. hyperfine states) 
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are split by an energy fkj a and coupled via optical transi- 
tions of frequencies uj c ±uj a /2 to an excited state |e) with 
population decay rate T (see Fig.[TJ:). We assume that a 
resonator mode of interest, with linewidth k, is tuned to 
the frequency uj c , i.e. it has equal and opposite detunings 
±A = ±uj a /2 from the transitions |f) -H> |e) , ||) -R> \e). 
For simplicity, we assume that the two transitions have 
equal strength and that all atoms are equally coupled 
to the resonator, with single-photon Rabi frequency 2g. 
The Hamiltonian is then given by 



H = flljJ r C 



tc+?l E(f [ItMtliH^aiJ+WcleMel, 

2 = 1 

+ 9[c\e) i {% + c\e) i {\\ i +H.c,\), (1) 

where c, are the photon annihilation and creation op- 
erators for the cavity mode and the index i labels the 
individual atoms. 

As we are interested in effects of the light shift and 
not in populating the excited state, we assume a large 
detuning |A| ^ n,T,g. For sufficiently low intracavity 
photon number (cJ c) <C (A/g) 2 , we can adiabatically 
eliminate the excited state [21| to arrive at an effective 
Hamiltonian for the dynamics in the two ground states: 



where S z = £ti(lt) 4 (t| 4 - (|| 4 )/2 and fi = 2 3 2 /|A|. 
(We shall later treat semiclassically the effect of spon- 
taneous emission from |e) into free space.) The first 
term in H e g describes the energy of the free field in the 
cavity, while the last term represents the energy of the 
atomic system with population difference 2S Z between 
the two states. The interaction term Mlc^cS z can be al- 
ternatively grouped with the first term to describe the 
shift of the cavity resonance by an amount QS Z due to 
the atomic-state-dependent index of refraction originat- 
ing from the transitions |t) <H> |e) and \\.) -O- |e), or with 
the last term to describe the light shifts zthilc'c/2 expe- 
rienced by the atoms in states |f) and |4-)- 

We adopt an interaction picture where the spin state 
vector evolves with the atomic Hamiltonian H a = 
hui a S z and the spin operator S evolves with the Hamil- 
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hflc'cSz + hui a S z 



(2) 



tonian Hq = h ( 
[c<t)(t),S(t)] =oh 
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Causality requires 



23] . It follows that S z is a constant 

of motion, while S+ = J2i=i It), (4-U = + which 
characterizes the phase on the Bloch sphere, evolves as 



dS'1 



ifrT 1 [flo,S+] = in$lJS\c 



(3) 



for any positive integer n. 

The equation of motion for the cavity annihilation op- 
erator c arises from the interaction both with the atomic 
system, through Ho, and with the bath modes outside 
the cavity, and can be derived with the standard input- 
output formalism for cavity fields (2^, [H[ : 



iuj r — iQS ; 
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where b\ , hi are the annihilation operators for the input 
fields from the two sides of the cavity (Fig. [TJ. The term 
iflS z describes the atomic tuning of the cavity frequency. 

We consider a monochromatic driving field, input from 
one side, in a coherent state \0) defined by b\(t) |/3) = 
-^/ne~ luJt (3 |j8) and 62^) = 0; where the frequency 
lo = uj c + k/2 is tuned to the slope of the resonator mode. 
In order to resolve the cavity resonance, we consider in- 
teraction times t 3> and we assume that the field 
is turned on adiabatically so that its transient behavior 
can be neglected. We consider an atomic system pre- 
pared initially in a CSS along the x axis satisfying 
S x (0) \ipo) = S \ipo)- To obtain a linear shearing of the 
CSS uncertainty region, we assume Ckso that 

the fluctuations of S z with variance AS Z = S/2 induce 
proportional fluctuations in intracavity power. 

The steady-state solution to Eq. 2] yields c(t) = 

e~ lcJt K/3/(v / 2(7— iu})) \/3) u , where the complex frequency 
7 = k/2 + i(oj c + nS z ) accounts for the leakage through 
the cavity mirrors as well as the phase evolution of the 
field. We substitute this result into Eq. [3] to obtain 
the time evolution of arbitrary powers of S+ in the spin 
subsystem, 



d(sr, 



dt 



if n {S z )(Sl)p, 
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where (G)p = (/3\ O \/3) denotes a trace over the field 
state. To lowest order in (0/«) \S Z \ < (Q/k)t/S/2 < 1 
and for n < y/S, 

f n (S z ) = nQ |/?| 2 (1 + n(i - l)n/« + 2(n/«)5»). (6) 

For easier visualization, after application of the squeez- 
ing light the spin is quickly rotated back about z, S + = 



S+(t)e 



-i/l(0)t 



by the angle fi(0)t corresponding to the 



light level for S z = 0. It is useful to introduce the dimen- 
sionless shearing strength Q = Spo(2fl/ k) 2 , proportional 
to the spin S, average photon number po = \/3\ 2 Kt/2 
transmitted for S z — in the time t, and square of 
the differential atomic phase shift per transmitted pho- 
ton 20//t <C 1. Equations [SHU then yield, in terms of the 
initial spin operators /S+(0) and S z , 

(S+)p = e^/ s S + (0) (7) 
(S 2 + ) fj = e -(i+i)Q/s e 2iQS„/s S 2 + ^ (g) 

Using Eqs. [THE] and Ref. we can evaluate various 
spin expectation values of interest for the atomic input 
state |<0o)- In particular, 

H 2 = Y + f-(Y-f)e-^s(Q), (9) 
(S V S Z + S z S y ) = (2S 2 - S) sin (J^j Gs(Q/2), (10) 

and AS 2 = S/2, where G s (u) = cos 25 " 1 ^/ S). The 
correlation term {S y S z + S z S y ), which is zero for a CSS 
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(Q=0), displays the quantum correlations between S y 
and S z induced by the ^-dependent spin precession. In 
the limit of a large ensemble, S > 1, and small phase 
\QS Z /S\ ~ Q/VS < 1, we have 

A5 2 «|(l + Q + Q 2 ). (11) 

The three terms in this expression represent the vari- 
ance 5/2 of the initial CSS, the variance increase due 
to photon shot noise QS/2 oc po oc t, and the vari- 
ance increase due to cavity-mediated coherent feedback 
Q 2 S/2 oc Pq <xt 2 . Thus, for Q > 1 feedback stretches the 
uncertainty region more quickly than does photon shot 
noise. This allows spin squeezing by one-axis twisting in 
the open quantum system, even though tracing over the 
output light field results in dissipative dynamics of the 
spin subsystem. 

To verify that the spin state is squeezed, we calculate 
the spin variance along the z axis after rotation about 
the x axis by an angle —a: 

AS 2 a = \ (v+ - y/vF+W* cos [2(a - a )]\ , (12) 

where V± = AS 2 ± AS 2 , W = (S y S z + S z S y ), and 
tan 2a = W/V- are specified by Eqs. IMTUl 

For moderate squeezing 1 <C Q 2 <C 5, the minimum 
and maximum variances, normalized to the CSS vari- 
ance, are ct 2 o « 1/Q and cr 2 Q+7r/2 Q 2 , where ct 2 = 

AS 2 / (5/2). While the antisqueezing increases as Q 2 due 
to cavity feedback (Eq. Hip , the minimum variance de- 
creases only as Q _1 because of the uncertainty in spin 
precession angle resulting from photon shot noise. The 
growth in the uncertainty product a ao a ao+7T / 2 ~ \[Q 
results from ignoring information in the outgoing light, 
which is entangled with the ensemble spin. 

For a given atom number 25, the squeezing improves 
with photon number po until the curvature of the Bloch 
sphere leads to reduced correlation between S y and S z 
[3, EH ■ Expanding A5 2 to lowest order in the charac- 
teristic phase variance Q 2 /(2S), we find <t 2 o = Q^ 1 + 
Q 4 /(2AS 2 ). The curvature of the Bloch sphere thus lim- 
its the minimum uncertainty to er 2 urv = (5/4)6 _1 / 5 5 _2 ^ 5 , 
reached at a shearing strength Qcurv = 6 1 / 5 5 2 / 5 . This 
is the same scaling with atom number as derived for a 
related scheme in free space (l8j . 

The attainable squeezing is also limited by photon 
emission into modes other than the mode of interest 
[I7L [25| . Squeezing requires two-fold light-atom in- 
teraction, namely the tuning of the cavity by the atoms 
in combination with the light shift by the modified intra- 
cavity intensity. As each process is proportional to the 
real part Re(a) oc f2 of the atomic polarizability a, the 
shearing strength Q oc O 2 oc Re(a) 2 oc A -2 oc Im(a) is 
proportional to the average number of photons 2r emit- 
ted into free space per atom during the squeezing process. 
Therefore, free-space scattering cannot be ignored at any 
light-atom detuning A, as is also evident if we express 



the shearing strength Q in terms of the single-atom co- 
operativity 77 = Ag 2 /(nT) as Q — ASr/r. This expression 
furthermore shows that achieving Q > 1 at low photon 
scattering probability, r -C 1 — as is necessary to main- 
tain coherence in the system — requires a large collective 
cooperativity, 5 77 ^> 1. 

For the symmetric level scheme considered here (Fig. 
QJ, Rayleigh scattering (where the atom returns to the 
same ground state) occurs at equal rates for states |f) , ||) 
and causes no decoherence since it provides no informa- 
tion about the atomic spin [26]. The random changes 
in S z arising from Raman scattering |f) — > \e) — > \l) or 
||) — > |e) — » |t), on the other hand, reduce the correlation 
between the time average S z during the squeezing, which 
determines the evolution of S y , and the final value S z (t) 
when the rotation about the a;-axis is performed. Replac- 
ing S z by S z in Eqs. [7] and [8] and evaluating the modified 
variance A5 2 and covariance (S y S z (t) + S z (t)S y ) in the 

large-5 limit using 2(S z (t 1 )S z (t 2 )) / S = e~ 2r ^- t2 \/\ we 
apply Eq. [T^] to calculate the minimum normalized vari- 
ance er 2 Q r . To lowest order in the number r « 1 of 
Raman-scattered photons, and ignoring curvature effects 
for the moment, 

As a function of r or time, the variance first decreases 
below the CSS variance due to the coherent feedback 
but eventually rises again when the noise added by pho- 
ton scattering into free space becomes appreciable (Fig. 
[5]) . Consequently, there is an optimum shearing strength 
Qscatt = y/SSr], at which the dynamic squeezing reduces 
the spin variance by a factor l/cr 2 or = y/3Srj/2 be- 
low that of the CSS. (This is the same scaling as in 
measurement-induced squeezing (l7l . H3|.) Squeezing is 
thus possible even for a resonator that is weakly coupled 
to a single atom (77 <C 1), as long as the collective co- 
operativity is large (Srj 3> 1). The limit A5 2 urv is only 
reached for 5?7 5 > 1; this is readily satisfied in moder- 
ately coupled resonators (77 > 1), but in free-space-like 
[IH situations (77 <C 1), for any reasonable atom number 
25, the squeezing is limited by scattering long before the 
curvature of the Bloch sphere becomes significant. 

By projecting atoms into |f) and Raman scattering 
not only adds noise but also shortens the spin vector. 
Provided 5 77 S> 1, however, this effect is negligible at 
the optimum squeezing point r opt = y/Z/{l§Srj) <C 1. 
Note also that it may be possible to find other transition 
schemes where Raman scattering is suppressed [25j|. 

Fig. [5]shows the minimum variance <r 2 r as a function 
of shearing strength Q for various values of the coopera- 
tivity 77, calculated from the full expressions of Eqs. IMTUl 
modified as described above to account for Raman scat- 
tering. The parameters 77 = 0.1, 5 = 10 4 are similar to 
those used to achieve squeezing in Refs. (IH \2l\. 

The three approximations of low saturation of the 
optical transitions, adiabaticity of the input pulse on 
the scale of the cavity lifetime, and small tuning of 
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FIG. 2: (Color online) Minimum normalized variance <j\ a ^ r 
as a function of shearing strength Q for S = 10 4 and various 
single-atom cooperativities rj = 0.001, 0.01, 0.1, 1 (solid lines). 
The dashed line shows the limit o"c Urv due to the curvature of 
the Bloch sphere when free-space scattering is ignored. The 
dotted line shows the variance neglecting both free-space scat- 
tering and curvature, scaling as 1/Q for Q > 1. 

the cavity by the atomic quantum noise are all realis- 
tic. Since the excited-state population e — (c' f c)g 2 /A 2 



is inversely proportional to the interaction time t, with 
ent = (k/ g) 2 Q / (8S), the first two conditions e <C 1 and 
Kb 3> 1 are always consistent. For S — 10 4 , r\ = 0.1, 
k = 2tt x 1 MHz, g = 2tt X 0.4 MHz, and A/T = 500 
[l5L \2l\. requiring e.g. e < 10~ 5 also ensures that the 
optimum shearing parameter Q « 50 is only reached at 
t > 400/k. The dependence of intracavity power on S z 
stays well in the linear regime, as fl^/ S/2 = 7 x 10~ 3 k. 

The expressions derived here well describe recent re- 
sults using a dilute ensemble of 87 Rb atoms inside an op- 
tical resonator in the moderate-coupling regime (27) . Un- 
like measurement-induced squeezing (l4l . [l5| , the present 
method does not require detection of the squeezing light 
and produces unconditionally squeezed states. Such 
states may be used in an atomic clock (28l . |29| to re- 
duce the quantum projection noise of the readout. Our 
scheme is applicable to both microwave and optical clock 
transitions, as it can be generalized to any configuration 
with different light shifts for the two ground states. 
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